A LIDAR SYSTEM FOR REMOTE PROBING OF THE 
LOWER ATMOSPHERE 


Prepared By- 

Charles D. Craig and Robert Bartz* 


Lars E. Olsson** 
Principal Investigator 


"E. Wendell Hewson 
Chairman, Dept, of Atmospheric 



Department of Atmospheric Sciences 
Oregon State University 
Corvallis, Oregon 97331 




. {NASA— OR- 141013) A LIDAR SYSTEM t?n b ~ ~ 

'REMOTE PROBING OF THE LOWER ATMOSPHERE 
(Oregon State Univ.) 53 p HC $4*25 

CSCL 171 


G3/36 


N75- 13241 


Unclas 

03653 


November 1974 


* Present Affilliation , School of Oceanography, Oregon State University. 

#* Present Affilliation, Swedish. Meteorological and Hydrological Institute, 

Stockholm. 

NOTE: The original version of this report was prepared in December, 1972. 



TABLE OP CONTENTS 


LIST OF TABLES 

LIST OF FIGURES 

1. SUMMARY . 

2. SYSTEM TECHNICAL DESCRIPTION . . . 

3. TESTING OF THE SYSTEM 

b. CONCLUSIONS 

5. SUGGESTED EQUIPMENT MODIFICATIONS 

6. SUGGESTIONS FOR FUTURE RESEARCH . 

7. REFERENCES 

8. APPENDIXES . . 



i 


Page 

LIST OF' TABLES 

Table 1. Lidar System Characteristics . . , , 4 

Table 2. Mixing Heights Estimated by Lidar Compared , . . l6 

with Estimates Made by Other Techniques 



-ii- 

LIST OF FIGURES Page 

Figure 1. Lidar System Schematic 2 

I 

Figure 2. Lidar System .....; 3 

Figure 3. Azimuth and Elevation Positioner with "i" Beam Mounted . 6 

Figure 4. Optical Receiver , 6 

Figure 5* Holoheam Series 300 Pulsed Ruby Laser System ...... 10 

Figure 6. Holobeam Pockels Cell Q-Switch 11 

Figure 7- Date Display Electronics Schematic Diagram 12 

Figure 8. Data Reduction Scheme 14 

Figure 9- Typical Lidar Returns Showing Top of Haze Layer .... 15 

in Western Oregon, Range 375 m/ cm 

Figure IQ. Range Corrected Lidar Signal and Temperature Versus . . 10 

Altitude for 1600 PST, 30 August 1972. 



A Lidar System for the Remote Probing of the Lower Atmosphere 
1 . Summary 

The Department of Atmospheric Sciences at Oregon State University has 
developed and tested a lidar* system under NASA Grant No. NGR 30-002-01*8. 

The lidar system employs a Q-switched ruby laser that radiates an in- 
tense pulse of light into the atmosphere, thereby illuminating aerosol particles 
in the laser beam. Light backs cat tered from these particles is received by a 
cassegrain telescope that is equipped with a sensitive photomultiplier tube. 

The output signal of the photomultiplier tube is displayed with respect to 
elapsed time (A-scope display ) on an oscilloscope and photographed on polaroid 
film. The film is manually digitized and subsequently processed by computer to 
obtain normalized apparent backscatter coefficient profiles which can be in- - 

terpreted (with some assumptions ) as approximate profiles of relative aerosol 
c oneent rati on . 

The system is capable of measuring aerosol profiles to ranges greater 
than three kilometers. 

2. System Technical Description 

The lidar system is illustrated schematically in Figure 1 and pictorially 
in Figure 2. The lidar system characteristics are summarized in Table 1. The 
system consists of an azimuth and elevation positioner, optical mount, a ruby 
laser transmitter, an optical receiver, electronics consoles, and a data dis- 
play and recording system. 


* lij 5 ht detection and ranging 
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FIGURE 1. LIDAR SYSTEM SCHEMATIC 








TABLE 1: LIDAR SYSTEM CHARACTERISTICS 


Transmitter 
Laser type : 

Mode of Operation: 

' Wave Length : 

Energy/Pulse: 

Pulse Width: 

Peak Power: 

Beam Divergence: 

Pulse Repetition Rate: 
Cooling: 

Receiver 

Telescope: 

Collecting Area: 

Field of View: 

Spectral Band Width: 

Light Detector: 

Response: 

Response Time: 

Current Amplification: 
Quantum Efficiency: 

Enclosure and Optical Mount 
Enclosure: 


Optical Mount: 

Data Display 
Oscilloscope: 

Camera: 

Data Processing 
Digitizer: 


Holobeam, series 300, Ruby, Q-switched, 
model 330 

Pockels cell Q-switched 
6943A 

Variable 0.2 to 2.0 joules 
10-30 nanoseconds 
Variable 10 to 100 MW 
3.8 milliradians 
6 PPM 

air cooled water- circulating system, above 
ambient temperature regulator 

Questar 89 IW field model 
62.2 CM 2 
5 milliradians 

10 nm, optics tech interference filter 

Photomultiplier, R.C.A. , 7326 

S-20 

-9 

2.5 x 10 seconds 

1 x 10 6 

2.5% at 6943A 


Optical system is mounted to a sturdy 
alluminum "I'* beam and completely enclosed 
by dust covers. Electronics are contained 
in two consoles. 

0-90° elevation angle, manually controlled, 
0-360° azimuth angle, motor driven. 


Tektronix, type 556, dual channel, 50 MH 
bandwidth 1 

Tektronix, Type C-12, 3% x 4 h polaroid film 


Calma 303, manual curve tracing to magnetic 
tape 

CDC 3300, data out in normalized relative 
values of apparent backscatter coefficient 
vs., range. 


Computer 



Azimuth, sind Elevation Posit loner 

The optical system (transmitter and receiver) can be manually positioned 
to any elevation angle from 0° to 90 ° and rotated electrically in azimuth, from 
0° to 360° with the positioner shown in Figure 3. Provisions are made to 
manually read or electrically transmit position information. The system is 
leveled by viewing a bull’s eye level attached to the positioner while adjusting 
the lock down pads. 

Optical Mount: 

A sturdy lt I M Beam is used to mount the laser transmitter and optical 
receiver. The optical system is enclosed by dust covers , with two ports on 
the front of these covers to enable transmission and reception of light. For 
safety purposes these ports are covered except during system use, preventing 
unintentional transmission of laser energy. A port is also provided on the 
receiver dust cover to view the finder scope used to observe the target area 
and to facilitate receiver alignment. The receiver axis can be moved ±5° with 
respect to the I-beam axis to accomplish alighment of the transmitter and 
receiver beams. 

The enclosure was designed to allow the transmitter beam to be transmitted 
coaxially. This modification will allow studies to be conducted at ranges less 
than 100 meters. 

Optical Receiver 

The receiver , which is shown in Figure k , was constructed by interfacing 
a Questar field model telescope (8p mm. aperature) and a Pacific Photometric 
photomultiplier assembly, model 62 F, with an adaptor designed by- Robert Bartz. 
The adaptor houses the field stop and narrow band interference filter. The 
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FIGURE 3- AZIMUTH AND ELEVATION POSITIONER 
WITH "I" BEAM MOUNTED 



FIGURE U. OPTICAL RECEIVER 
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receiver is mounted on a base assembly which, permits azimuth and elevation ad- 
justments required for alignment of the receiver and transmitter optical axes. 
A reticle in the telescope eyepiece is used to align the optical system and to 
aim it accurately at a known target. The telescope is also equipped with a 
wide angle finder scope that permits safe operation of the system by enabling 
the operator to view any obstruction that may come between the system and the 
intended target. Technical data for the receiver is presented in appendix A. 

Laser Transmitter 

The transmitter system components, shown in Figure 5 and 6 include a 
Holobeam Series 300 pulsed ruby laser system (laser head, end mirrors, optical 
rail, remote station, power supply) a diy Pochels cell Q-svitch assembly, and 
shutter electronics. The Q-switched laser emits radiation at 0.6^3 microns 
with power levels exceeding 100. megawatts. The output power of the laser is 
monitored through the rear mirror of the laser with a modified Quantronix 

model 506 power monitor and displayed on one channel of a Tektronix model 556 
dual channel oscilloscope. 

Electronics 

H>e laser is fired by closing the flash contacts on the oscilloscope 
camera at a repetition rate determined by the operator. The laser power supply 
and its safety circuits are contained in the power supply console cabinet. The 
console operates from a 120 volt 60Hs source. The remote station is connected 
by a 10-foot multiconductor cable to the power supply. The voltage to which 
the power supply capacitor bank is charged determines the laser output energy 
and is preset at the remote station. The energy required to fire the flash 
lamps to start laser pumping is initially stored in the capacitor hank and then 
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released through, a aeries inductance upon command from the firing circuit. 

The laser cooling unit is contained within the power supply cabinet and is 
accessible through the rear panel door. 

The laser Q— switch (Figure 6) consists of the Pockels cell assembly, the 
shutter electronics power unit, remote control station, and connecting cables. 
The Pockels cell assembly contains a dry KD*P crystal, a polarizer, and a maxi- 
mum reflectivity dielectric mirror. At the instant the flashlsmp in the laser 
head is fired, this cell is cross-polarized with respect to the polarized 
light from the laser rod. In this condition, the cell is off and lasing action 
is inhibited. After a set delay time (controlled by the Pockels cell remote 
station), during which the population inversion of excited atoms within the 
rod rises to extremely high levels, a high-voltage pulse is applied to the 
crystal. This causes a 90 degree rotation of the plane of polarization of the 
cell and lasing occurs as a giant pulse of 15-30 nanosecond duration with peak 

g 

power exceeding 10 watts. Before conventional lasing can begin or a second 
pulse can occur, the voltage applied to the Pockels cell crystal decays and the 
plane of polarization of the crystal returns to its blocking condition. 

Optical Receiver 

The field of view of the telescope is limited to five milliradians by a 
field stop. After passing through the field stop, light collected by the 

o 

telescope passes through an interference filter with a bandwidth of 10A to an 
RCA 7326 photomultiplier tube. Photomultiplier output is buffered by an 
emitter follower and applied to one channel of the oscilloscope. 

Receiver gain is controlled by adjusting the photomultiplier tube supply 
voltage. Sensitivity and current amplification characteristics of the photo- 
multiplier tube are given in appendix A. 
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The photomultiplier power supply is a Pacific Photometries Model 203 
high voltage power supply. Power supply voltage can be adjusted from zero to 
-2000 volts D.C. This power supply also contains a ±15 V.D.C. output that is 
used to provide voltage to the emitter follower used to buffer the photo- 
iffultiplier tube output* 

Display and Recording 

The laser output power and atmospheric backscatter return signal are 
displayed on separate channels of a Tektronix 556 dual channel oscilloscope. 
These signals are recorded by photographing the oscilloscope display with a 
Tektronix C12 oscilloscope camera using Polaroid fil m . 

The camera shutter flash contacts, which are controlled remotely by the 
observer operator, ar*e used to fire the laser system. AH controls can be 
preset to display, record and fire the system which allows one man to operate 

the lidar system. The data display electronics are shown schematically in 
Figure 7. 

Electronic Calibration 

The lidar system is calibrated in range only by means of a Tektronix 
l8l time mark, generator. 

Data Reduction 

fte retUrn Signal is displayed on an oscilloscope and photographed with 
a Polaroid camera. Return signal strength and range Information will be ex- 
tracted from the photograph as x, y coordinates by means of a Calma 303 
digitizer (available at OSU Computer Center) which stores the coordinates 
directly on magnetic tape. The x, y coordinates of the photograph and the 




FIGURE 5. HOLOBEAM SERIES 300 PULSED RUBY LASER SYSTEM 
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required system parameters will then be fed into the OSU GDC 3300 . computer 
The computer output will consist of range-corrected, normalised, relative 
values of the apparent backscatter coefficient vs. range. A block diagram 
of the data reduction scheme is shown in Figure 8. 


3* Testing of the System 

During August and Septunber of 1972, . a special field study to investigate 
the effect of agricultural field burning on the surface energy budget was under- 
n near Corvallis, Oregon. During the period from August 28, 1972 to 
September 1, 1972, several vertical aerosol profiles were obtained by operating 
the Ildar from a site on the Oregon State University Campus. The vertical 
aerosol profiles obtained with the lidar system were analyzed and the "actual 
mixing heights" (defined here as the height at which there is a pronounced 
decrease in aerosol concentration 1 ] were compared with mixing heights 
estimated frcm data collected by other means during the above mentioned 

P ‘al field study. Figure 9 shows typical polaroid photographs of the lidar 
return signal obtained during this period. 


Mixing heights obtained fran the various observational methods are pre- 
sented in Table 2. It should be kept in mind that the radiosonde observations 
were made in Salem, approximately 30 miles north of Corvallis, while the air- 
craft ascents, pibals winds aloft measurements, and acoustic soundings were 
made approximately 5 miles south of Corvallis. There is a general agreement 
between the various simultaneous estimates of mixing height on most occasions. 
Nevertheless, there are considerable spatial and temporal variations in the 

estimated mixing heights, a finding in agreement with the results of previous 
investigations. 1,2 ’ 3 
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FIGURE 8. DATA REDUCTION SCHEME 
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FIGURE 9- TYPICAL LIDAR RETURNS SHOWING TOP OF HAZE LAYER IN 
WESTERN OREGON, RANGE 375 m/cm 
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Mixing Heights (Meters ) Estimated by 


Sounding Time 
(LST) 

TABLE. 2 

Lidar Compared with 
Other Techniques 

Estimates 

Made by 

L-R 

RAOB AIRCRAFT 

w-s 

A-S 

8-28 

1645 

1200 

1400 

1400 

— 


1945 

1200 

— 

— 



2010 

800, 1200 

— 

— 

200 


2400 

600, 1100 

— 

1100 

200 

8-29 

0400 

270, 600 

740 


100 


0630 

340, 890 

— 




0730 

420, 890 

— 




1115 

495, 940 

— 




1500 

730 

1700 


700 


1715 

600 

— 

900 

250 


1945 

370, 670 

— 


400 


2200 

100, 340, 700 

— 

500 

450 


2400 

330 

— 


450 

8-30 

0400 

850 (cloud) 

100, 750 

1500 

200 




2650 




0700 

340, (670) 

— 


100, 300 


1000 

700 

— 700 

400, 750 

700 


1400 

800, 1000 

— 

990 

400* 


1600 

890 

1250, 2100 

1200 

400* 


1945 

890 

— 


400* 


2015 

270, 670 

— 


400* 

8-31 

0400 

290, 440, 730 

100, 530 


150 


0700 

760 


750 

300 


1000 

600 


800,-900 

600 


1330 

800 


1000 

400* 


1600 

760 

1460 

1100 

— 


1900 

320, 440, 650, 690 






2200 

310, 730 


800 

450 


2400 

490 



500 

9-1 

0400 

130, 180, 240 

110, 370, 500 




0800 

440, 730 



• 


*) Where no clear indication of stable layer aloft was seen on the A-S chart the 
estimated value represents the height of the .tops of the thermal plumes. 

LR — lidar 

RAOB — U.S. National Weather Service Radiosonde 
A-g- Sounder 

W-S- Win as hear Layer from Pibal Measurement. 
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Figure 10 shows the range-corrected lidar signal and the vertical 
temperature sounding obtained at l600 PST on August 30, 1972. There is a pro- 
nounced decrease in range corrected signal return through a layer near the 
base of an elevated inversion. 

Conclusions 

The results presented in the preceding section indicate that the opti- 
cal radar system developed at OSU under NASA Grant No. • NGR-38-992-048 meets 
the design objectives set forth in the Development Plan. 1 * Accordingly, the 
instrument is capable of gathering data on a spatial and temporal resolution 
that is impossible or infeasible with most other techniques. This capability 
enables one to construct a detailed picture of atmospheric structure which may 
furnish valuable insight into the nature of some atmospheric porcesses. 

5 * Suggested Equipment Modifications 

The following modifications and additions to the instrument system are 
suggested: 

1. The data display system should be converted from oscilloscope Polaroid 
photography to a real time display that presents in raster or map form the 
spatial and/or temporal distribution of backscatter profiles. This would 
make possible RHI, PPI and time-height cross section type displays. 

2. The data recording system should be converted from oscilloscope polaroid 
photography to a fast analog to digital converter that samples and digitizes 
the return signal. Integrated circuits that are fast enough for this purpose 
are commercially available. Once the signal has been sampled and held in digi- 
tal form, it could be immediately analyzed with a small on-line computer and 
stored on a magnetic disc recorder for instant replay in the display type of 
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FIGURE 10- RANGE CORRECTED LIDAR SIGNAL AND TEMPERATURI 
VERSUS ALTITUDE FOR 1600 PST, 30 AUGUST 1972 


Temperature (30 -August 1972) Salem RAOB 1600 PST. 
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the operator's choice. 

3. The system should he made mobile by mounting it in a small truck equipped 
with its own power supply and air conditioning equipment. 

The transmitter and receiver optics should be made coaxial to permit 
measurement of signal returns from target volumes located at distances less 
than about 150 meters from the instrument. 

These modifications would significantly enhance the present capability 
by permitting mobile operation, observation of the close-in return, on the 
spot inspection of the atmospheric area being surveyed, and real time data 
analysis. It has not been the purpose of this section to describe in detail 
the electronic and mechanical constructs necessary for these modifications. 
Allen and Evans 5 have reported details of a recently developed system that 
in co npo rates most, of* these modifications » 

6. Suggestions for Future Research 

Although the development of an active observational research program 
involving the optical radar system has been delayed by lack of funds, the 
following uses for the instrument have been seriously considered: 
l) Observations of the spatial and temporal variations of the mixing 
height related to the influence of topographical inhomogenities on 
wind fields and surface heating, dynamic and thermodynamic influences 
of local wind systems, radiative effects of polluted layers, and urban 
heat island effects. 

2j Observations of the changes in the stratification and transport of sus- 
pended particulate matter during the development of mesoscale thermal- 
tidal wind systems • 
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3) Observations of plume behavior and diffusion in mountainous terrain 
with particular emphasis on the behavior of a plume approaching an 
obstacle such as a butte or a bluff. 


4) Observations of sea-salt spray generation in the coastal surf zone. 
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APPENDIX A 


• OPTICAL RECEIVER TECHNICAL DATA 



A-l 


QUESTAR SPECIFICATIONS 


TYPE: Maksutov Cassegrain Catadioptric. No coma, astig- 
matism or spherical aberrations. 

CLEAR APERTURE, 3.5 inches, 89 mm. 

FOCAL LENGTH, Basic Visual, 50.5 inches, f/14.4, 1300 mm. 
FOCAL LENGTH, Camera Close, 56 inches, f/16, 1400 mm. 
FOCAL LENGTH, Camera with Ext. Tubes, 64 inches, f/18, 
1600 mm. 

FINDER LENS, 4' FI., 4x and fix. Field 12° and 8°. 
POWERS, 40~80x eyepiece, Field 55 min. 

POWERS, 80-160x eyepiece, field 42 min. 

POWERS LIMIT: Resolves 1 sec. arc at 50 feet EFL. 

FIELD OF VIEW, photographic model, t°30 min. 


EYEPIECES, 26mm. Koenig, 50°Ap. Field; 12 men. 5-lem 
Erfle, 75°Ap. Field. 

AMPLIFYING OR BARLOW LENS, minus 43.9 mm. FL, 

ERECTING PRISM, Star Diagonal rype, 90®. 

STANDARD SUN FILTER, 1-1/2' aperture glass fiat, off-axis 
mount with exact measured thickness of chromium coating 
allowing a .002% to .003% light transmission. 

FULL-APERTURE SUN FILTER, 3-1/2' clear aperture, ex- 
act measured thickness of chromium coating allowing a 
.002% to .003% light transmission. 
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APPENDIX B 


LASER SYSTEM TECHNICAL DATA 
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HOLOBEAM, INC. 


TECHNICAL DATA SHEET 


SERIES W 1 ™ POCKELS CELL Q-SWITCH ACCESSORY 

CUSTOMER SEfrFfe HOLOBEAM JOB NO. VSAOOfr 


SHIP BY 


LASER SYSTEM DATA 

» 

LASER ROD 
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OUTPUT REFLECTOR 
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APPENDIX C 


THE LI DAP RANGE EQUATION AND COMPUTER PROGRAMS FOR ITS SOLUTION 


The following text and computer programs were supplied by Dr. G.W. 
Betlike, Space Sciences Lab., General Electric Company, P.0. Box 8555, 
Philadelphia, PA 19101, and are reproduced verbatim. 



APPENDIX A 

THE LIDAR RANGE EQUATION AND COMPUTER 
PROGRAMS FOR ITS SOLUTION 

^ • Range Equation Derivation 

It is assumed that the laser light source is pulsed,with light 
pulses which are short in length compared to the measurement ranges 
of interest. It is also assumed that the laser beam is collimated to a 
divergence which is smaller than the lidar receiving system angle of 
acceptance (field of view). It is further assumed that the transmitted 
laser beam is coaxial or nearly coaxial with the receiving system, such 
that the propagating laser beam is located entirely within the receiving 
system field of view during the period (ranges) of interest. Because the 
laser beam is located entirely within the receiver field of view, only 

the total power (P L ) and energy (E) of the laser beam need be considered, 
and not the beam intensity* 

The total power of incident' light (P.,) scattered in direction 8 by 
illuminated scattering length {&) of gas is 

P s (e) = P i ^k s { 0 ) (!) 

where k g (0) is the directional scattering coefficient. * 

Now referring to Figure A-l, we see that a pulse of light which 
starts from the laser at time t Q and continues to be emitted until time 
has a length L, where 



RANGE (r) 


C-2 



TIME (t) 


Figure A-l. Laser pulse range vs. time diagram. This diagram 
illustrates the length (L/2) o£ the lidar system laser 
pulse scattering volume which is observed at time t^» 
for a laser pulse of length L. 
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L - c (t L - t 0 ) = c t L (2) 

As this short light beam of length L propagates away, light is continuously 
scattered back to the lidar system. However, Figure A-l illustrates 
that the scattered light observed by the lidar system at any one instant 
in time, t r , is scattered only from incident light located between ranges 

r l an< ^ r 2 which have a distance difference of L./2. Thus the illuminated 
length of gas that scatters light back towards the lidar system at any 
one time, t r , is 

^ = r 2 ' T i = L /2. (3) 

The light which left the lidar laser with power arrives at the 
range of interest (r) with reduced incident power due to atmospheric 
extinction. Thus 

P i = P l ex P t- rk e } ( 4 ) 

where k g is the total extinction coefficient of air from all sources 
assuming extinction to be constant with range. Finally, 

P L = E=E/t L {5) 

where E is the total laser energy per pulse. Now combining equations 
I through 5, we have 

P s (9) = (cE ■ k s (9)/2) exp (-rk e ). (6) 

The power of back-scattered light collected (P c ) by the lidar 
receiving objective is 

P c = P s ( 180 °) • T exp (-rk e ) 


(7) 



where Tis the solid angle subtended by the lidar receiving objective from 
the scattering range r. But 

r = 27T{1 - cosy) «jry 2 « A/r 2 (8) 

where Y is the angle subtended by the radius of the lidar receiving objective 
from scattering distance r, and A is the area of the receiving objective. 
Since it has been assumed that the light pulse is short in length compared 
to the measurement ranges of interest, we can assume with negligible 
error that F is the same at any one instant (t r ) from all portions of the 
illuminated scattering volume. Finally, the power of scattered light 
which reaches the detector (P^is 

P d = E o P c (9) 

where E Q is the total efficiency of all optics in both the transmitting and 
receiving sections of the lidar system* Thus combining equations 6 through 
9 we have 

P d = f E-k s (180°).exp (-2rk e ) . (10) 

If is in watts and the detector has a sensitivity S in amps /watt, 
then the voltage drop, V, produced by detector output current I across 
detector load resistor R is 

V = IR = P d SR (11) 

The directional scattering coefficient consists of both molecular 
(Rayleigh) scattering, and particulate or aerosol (M.ie) scattering. 

k g (0) = n(r) • otr (9) + C^mtr) • cr M (0) (12) 



Here, n is the molecular concentration (density) in air, <7^(0) is the 
directional Rayleigh cross section, m is the particulate concentration in 
a * r ’ i s the suitably size-integrated directional Mie theory-calculated 

cross section for particulates in the air, and is a factor which corrects 
the Mie theory-calculated <7 j^(9) to the correct value for the real atmospheric 
aerosol. Similarly, the total extinction coefficient is defined as molecular 
(Rayleigh) and particulate or aerosol (Mie) total scattering plus any 

s 

molecular absorption at the specific wavelength of interest. 

k e = n(r) • aR + C 2 m(r) - a M + k a (13) 

Here a R is the total Rayleigh cross section, kg. is the absorption coefficient 

of air at the laser wavelength, is the total Mie theory-calculated 

cross section, and C 2 is a factor which corrects the Mie theory-calculated 

<7 m t0 correct value for the atmospheric aerosol. As is clear from 

equation 13, k g is not constant with range. Consequently, we must 

replace the expression (rk g ) in equations 4, 6, 7, and 10 with 

i* r 

r k e -*a R j* n ( r )dr + C^ cr M J m(r)dr+rk & . (14) 

° o 

Finally, combining equations 10, 11, 12, and 14, one obtains the 
range equation 

V = cR SE Q — j j^n(r) • a R (180 ) + C ^ m(r) • a ^ ( 180 °)j 
exp - 2 Jcr R J^n(r)d r + C 2 7^J m(r)drj . (15) 

Here we have dropped the absorption coefficient of equation 14 because 
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we assume that the laser wavelength does not coincide with a significant 
atmospheric line or continuum. 

The atmospheric density, n can be expressed analytically if one 
curve-fits to a suitable standard atmosphere. For lidar system elevations 
that do not cover a range of more than a few thousand feet above sea level, 
the following expression is useful to about 50,000 feet scattering altitude: 

n = n 0 • |exp (-2. 3026 sin <p) -Z" r^sin^<pj {16) 

Here, n Q is the value of n at temperature T q and pressure p Q , while T 
is air temperature and p is air pressure near the lidar system, (p is the 
lidar system angle of elevation above horizontal, and Z' and Z" are 
suitable constants. If the lidar system is at sea level and r is in cm, then 
setting Z' = 2.31 x 10^ cm and Z" = 3.0 x 10 ^ cm ^ gives n values 
which agree with the AFCRL 1959 standard atmosphere to within 50.6% 
and <3% at 530,000 feet and s50,000 feet altitudes above the lidar, 
respectively. If the lidar system is at 2000 feet altitude, then the same 
assumptions give agreement to within ^ 1% and <6% at ^20, 000 feet and 
* 50,000 feet altitudes above the lidar, respectively. Of course for 0 << 

90 , one is not concerned with such high altitudes even when r is large. 
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II . Computer Solution of the Range Equation 

As has been pointed out in section I of this appendix, for the general 
case of non,- constant m(r), the range equation is an integral equation 
requiring a numerical rather than analytical solution technique. For 


purposes of clarity in this discussion, the range equation (see equation 15) 
is reproduced here in a form and notation which are consistent with that 
used in writing the computer program. It is assumed that Cj = C 2 = 


C = C\ 


where 


and 


with 


V = "T“ [ c '* ^( 180 °) • m + a ( 180 °) • n] 

r ^ 

/ r 

•exp(-2^J o n d r-2 c '^; om dr) 

K = (1/2)R SE q cAE 
n = x i [ x 2 r2 + ex P ^ x 3 r 0 


X 1 = n o T 0 P/ (T P o ) :>0 
X 2 = -Z" sin^ <p < 0 
x 3 = -(2.3026 sin <P)/Z' 

Performing the integration of n with respect to r, one can rearrange 


quation (17) into the following form: 

r 

y i ex P (y 2 ) = (m + y 3 ) exp (y 4 J mdr) 


(18) 



where 


Yj s Vr^/ [KC / c^ I ( 180°) ] 

Y 2 ~ 2ct r x i ^ ( x 2 r3 ^ 3 ^ + exp (x^r - 1)] 

y ? = <t r (I 80°> • n/ [C '• a M (180°)] 

y 4'- 2C ' a M • 

In practice the data points used for the computer solution of equation (17) 

are V.^ and r^ 1 =i = 200 with no restriction on the interval size other 

than the physically sensible one of taking fine intervals where V varies 

rapidly. Here, 2 0, and is unknown and not necessary for solution. 

The method of solution for equation (18) assumes m = m = 

o 1 

constant for the range r to r,. This then gives 
r^ o I & 

4 mdr = m i -< r i -%> * s j 

and allows equation (18) which is trancendental in m to be solved iteratively 
for m by Newton's method, reference (1). Two values of ix^ will satisfy this 
equation for a large range of values of the physical quantities involved. 

Only one of the values ofm 1 will be correct and will lead to the correct 
values for m. for all i. The existence of these two solutions can be 
explained best from a physical viewpoint. If one considers receiving a 
backscattered lidar signal from close range (rj), its level can be 
relatively small for one of two possible reasons? (a) either m^ is small, 
or (b) m^ through m^ is large, resulting in significant attenuation for 
range r^. The first solution (smaller rrij ) has been found to be the 
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correct one for £ 10 (igm/m with r^ - 4 x 10 cm, both of these 

values being much larger than normally encountered in use of the lidar 
system. Therefore, unless the physical situation (such as dense fog) 
dictates otherwise, the first solution is selected, although the second 
solution (if it exists) can be examined readily by using a simple program 
instruction. 

The remainder of the are obtained by successively using the 

previously calculated values of m. to evaluate the exponential term 

in equation (18) in the following manner: 

J r i mdr = S i . !+ [(m. _ j + m.)/2]- [r. - r. J 
o 

-S. (i>l) (19) 

Physically, this has the effect of using the average value of m(r) over 

the interval in computing this term. Combining equations (18) and (19), 

one obtains the final form of the range equation as treated by the compute 

program in solving iteratively at each rz (i > 1) for m using Newton’s 

i 

method . 

( m i + Y 3* CXP ^ y 4 * r i " r i- 1* m i ^ 

- Yj ex P fy 2 - Y 4 [S._| + (r. - r i . 1 ) m i _i / 2 ^ < 2 °) 

The computer program for performing this solution to the range 
equation is called RANGM and has been written in Fortran for desk side 
terminal service to the General Electric 605 computer. The RANGM 
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program requires the following inputs with the units indicated in parentheses: 

R (ohms), SE q (amps /watt), A (cm 2 ), E (joules), ^(180°) (cm 2 /(ngm sr)), 

T ( R), p (inches Hg), Z' (cm), <p(degrees), Z"(cm -2 ), C (unitless), o (cm 2 /pgm), 
number of data points^and the data input r (cm) and V (volts). For each 
input point, the program then prints out r (cm), V (volts), and m ( Mgm /cc). 

A complete listing of this program is shown below. 



READY 

SLIST RANGM 


RANGM 



s 

o 


u 

ur 


©e 





O 


V 


UJ 


fid 




01/13/70 17.373 


1 3 

333 20 
'■>333 33 
333 3 5 
33343 
330 53 

m 

330 73 

330 80 
03390 
33130 
3a 1 13 
33123 

33 1 33 
001.43 

30 1 50 
031*55 
33 1 63 
33 1 73 
33 1 83 

33 1 93 
30230 

332 1 3 
33220 
33233 
33240 
332 53 

302 63 
03261 
302 73 


D I M EM S 10 N VC 233 ) > RC203 ) 

REAL" 'J0*M,*,M'p,c1P> 

A SC 1 1 . M S (? 1 * M 5 G2 j 3 S G3 > M 5 6 4> F 1 L E > F I R S T 
CALL SET ERR ( 1 5) 

3SG1 = " ERR 1 " 

MSG2=" ERR2" 

MSG3=” ERR3' 1 ' 

M S 64= " ER R4 '* 

C=2» 99 71 E+ 13 
SGR= 1 . 6692 E- 27 
SGRT=1 .9925E-28 
N3 = 2 • 6871 E+ 1 9 
13=491 .7 


P3=29. 921 

PRINT: " R , S E3 » A * E^ SGM j S GM T" i R EA 0 : C APR » S E3 < A * t* SGM » S GMT 
PRINT: ”T>P*ZP* PHI .* ZPP> CP”i READ : T» P-> ZP* PHI* ZPP» CP 
- 13 PRINT: "SHI OLD 1ST* SLN DOM I NAT E? "i R EAD : FI R$T 
PRINTS "IS. DATA SOURCE ON. FILE?”: READ: FILE"" 

I FC FIL-E. E3 ."NO" ) GO TO 11 
REA0C"FILE1"> 1 ) NT 

1 FORMAT < I 5> 

R EA P < " FI L E 1" >2) <R( I )»V(.I)» I = 1 > N T ) 

2 FORMAT C El 6* 8* El 1 • 3 )’ 


1 1 


l 2 


GO TO 12 

PRINT : "NO • DATA P TS • " » R EA D : N T 
PR I NT: "RANGE AND VOLTAGE" 
PRINT:" " 


READ: < R C I ) j VCI >, 1= WNT> 
PRINT if” " 

CONTINUE 


EE PRODUCIBILmO|jHE 
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o 


u 

Ul 









*3329*3. 
30333 
3*T3 1 3 
3331 I 
03 3 1 2 
33323 
33333 
33 343 
33353 
33363 
00373 
30380 
30 393 
03 43 3 
00 410 
30423 
30 433 
30 4 40 
3*0 4 50 
00 463 
03.473 
33 483 
00490 
33 503 
03 510 
30 523 
00 530 
03 543 
03 550 
00563 
03 570 
30 58*0 
0*3 593 
03630 
00 6 1 0 
33623 
30 630 
00 640. 
30 65*0 
00 663 
30 670 
03680 
03 685 
■03 69'0 
03700 
m 7 1 0 
00723 
30 743 
037 50 
00 7 63 


K=CAPR/2. *SE3*C*A*E 
PHI=PHt/S?*2958 
SiMP=SI'J CP HI > 

TCST=Aas(SIMP3 

I FA TEST. LT* 4- * E- 1 0 > SINP-i * Et 1.0 
T1 =CP*SGMT 

T2=2 ♦ * < S GR*i>J 0 > * T0 /T*P/P3 
T3=SINP**2/3., 

T5=- 2 * 332 6* S INF 

74--ZP/75 

T6=< SGRT*N0 >*T0/T *P/P0. 
Y4=-2v*CP*SG0 
S=3 • 3 
-PRINT 3 

3 FO RM AT < / / 6 X, "R"> 11 X# ” V", .MX* ”M ” / / ) 
DO 1 30 1=1, NT 
Y 1 = ( V ( I > / T ! ) * ( R ( I > 

EX= T 5* C R ( I > / ZP > 

TFST=ASS(EX) 

.1 ,F (TEST • GT • ;*'0 1 > GO TO -40 

5110=0.3 

TRM=E* 


F I = 1 - 


3*3 SUM=SUM + TRM. 

TESTERS CTRM'SiM ) 

I F C TEST. L E. 1 . E-0 6) GO -TO 50 
FI=.FI + 1 * 

IF {FI- GT .20 • OOCALL' ERROR C 4,MSGl ) 

TRM=TRM/FI*EX. 

GO T O 3(3 

43 SIW= EXP < EX)- I * 

53 : Y2= T2* C ( - ZPP *R C l >*.* 3 ) * T3- T 4* SUM >' 

Y3=T6/T1 * C3-*T3*C-iPfi*R< I > **2 >+ EXP < T 5* < R < I >/ZP> > > 
I Ft I . EQ'. 1 ) DELR= R< 1 ) 

I F< 1 . GT. 1 > DEL R= <R ( I >-RU - 1)0/2* 

I F( I . E0. I ) C1=Y1 *EXP CY2) 

I F( I • GT* 10 S=S+**DELR 
I F< I * GT* 1 ) C 1 = Y 1 * EXPX Y2 - Y 4* S > 

C2= Y3* 


60 


C3=Y4*DELR 
MP = - (1 ./C3+C2) 


MPP=MP-1 ./Q3 
J=1 

IFCMP. GT*0* > GO TO 60 
IFCC1.GT.C2> CALL ERROR C 4, MSG2 ) 
M=3*‘ 

I F (MPP...GT * 3 * > M = MPP 
GO T) 70 

FMP= (MP+C2 >* EXP ( C3+MPO 
I FCC1 » GT • FMP ) CALL ERROR ( 4, MSG3 > 


44 ^3* 
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& 


C-12 



33773 

0.3 7 83 
33 790 
33830 
30 813 
33823 
00830 
33840 
30 8 S3 
00863 
33870 
03883 
30890 
00900 
33 9 1 3 
33923 
00933 
30 960 


IFCFTR5T. EQ. ,? N 0") G) TO 6 5 

M = 0. 

IFCCl.LE.C2) GO TO. 83 
GO TO 73 
65 M=Mpp ~ 

73 FM= CM + C2)*EXPCC3*MX. 

FpM= (1 .+C3*CM+fc2> )*EXPCC3*M> 
SAVMsM 

M=SAVM + CC1-FM j/.FPM 
TF.ST=A3S< CM 1 SAVM >/ CM+SAVM) ) 
J=J+1 

IFCJ.GT.50> CALL ERROR < 4»MSG4> 
1 FC TEST* GT'» 1 • E-0 6) GO TO 73 
80 S- S+M + D.ELR 


103 PR INT 4* R C I ) f V C I )-»M’ 
4 FORMAT C1P3E12* A) 

GO TO 10 


EN D 


In order to test the RANGM program listed above, it was found 
to be convenient to write an additional program called RANGV which will 
compute m. given V, and r.. This program simply solves equation 
(17) for V. using the procedures outlined in the previous discussion of 
RANGM for the evaluation of the integrals. This computation requires 
no iteration or choice of solutions since V\ is uniquely specified by 
m. and r^. The input is the same as for RANGM with the exception that 
the data input are now r^ and m . The program prints out r^, v\» 
and m. the same as RANGM. A complete listing of the program is 

shown below. 

READY 

*L I 5 l RANGV RANGV 

01/13/70 17.302 

•000 10 DIMENSION MC230 )> RC203 )> VC230) 

03020 REAL N0 * M> K 


013 





33331 
04343 
333 53 

3m 6 3 

333 73 

333 S3 
333.93 
'• 133 

m t i a 

my?, a 

33 T30 
33 1 43 
33 .2 1 3 
33323 
33221 
33243 
03253 
002 60 
33273 
00271 
302 72 


33283 



© 

u 

UJ 



33293 
00330 
00 31 0 
33 323 
33333 
33343 
30353 
333 63 
30 3 73 
03383 



■00390! 
00 433 
33 4 13 
03423 


© 

/ 

u 

UJ 

O' 


00433 
30 4 43 
004 50 
03 4 63 
03 4 70 
33480 
3,0490 
30495 
03 513 
03520 
30 530 



33 540 
33 5 53 
03 560 
•03 570 


ASCI I 05 81 1 RTFL 
CALL SETERRCT5) 

■MSG 1 ERR l‘* 

C=2- 99 71 E* 1 3_ 

SGR=1. 6692E-27 
SGRT= 1 » 992 5E-28 
N3 = 2.'6871 EM 9 
T3= 49 1 ♦ 7 ' 

P3-29* 921 ............ 

P2IMT: "R» SE3» A* E> SGrt* SGMT" i READ: CAP R» S£0> A# F» SGM* S.GMT 
PRINT: ”T» P* ZP» PHI '» ZPP* CP rt J READ: TiPvZp, PHI t ZPP* CP 
10 PRINT: "WRITE A 3CD FI L E?" » READ: RTFL 
PRINT:*’'! ). DATA pfs.’*J REAtJsNT 
PRINT: "RANGE AND MASS” 

PRINT:'* ” 

READ : C R C I ) » MCI) * I = l > NT* 

K=CAPR/2« *SE0*9*A*-E 

PHI=PH 1/57.2958 
S1MP=SIN(PHI ) 

TEST=A8S I SINP4 
I F (TEST »LT • 1 • EM 0 > 
n=CP*SGMT 
T2=2»* (SGR*N0) *T0>T 
T3=SINP**2/3«. 

T,S=^2.302 6* SIMP 
T4=- 1 ZP/T5 

T6= CSGRT *N3 > *TQ/T*P/P® 

PRINT 2 ' 

2 . FORMAT C//6X* ”R"» l 1 X* ” V*% I I X* "M”// ) 

D? i m - i = i*nt' 

IFCI. EQ* D S=M ( 1 ) *R,< 1 > 

I F C I • N E* 1 > S- S+ CMC I )’+M C-I - 1 ) J /2 * * t RC I > - R < I - f> ) 
EX=T5*CR< I >4ZP> 

TEST- ASS CE>0 

I FC TEST* GT~ .01 KG? T) 43 ' 

SIM-3. 0 • 

TRM= EX 
FI = 1 . 

30 SUM= SlWt-TRM 

T EST= A 8 S C TRM / S l H > 

I F C T EST- L E. 1 . E- 3 6 > G? TO 50 
FI = Fl+l-. 

fFCFI. GT.20. ) CALL ERR0RC4,MSGl > 

TR M= TRM/ F I * EX 
G? T ) 33 

40 SUM=EXPCEX>“t,« 

53 TRM=T2^2 • * C C-ZPP + R C I j **3 ) *T3-'T4* SUM )■ 

TRM = C TRM +.C ?* 3 GM * S 
TRM = EXP (r 2.* TRM) 

FCT=KXRC I ) **2* < T 1 *M C I ) +T6* C C-ZPP*R< l )**2)*3« *T3 
+ EXP C T5* C R C I )/ZP) ) ) > 




RECOMM 


C-l 1 * 


•M5R0 
m 59 5 

1 4 3 

V< I > = FCT*TRM 

PRI 'JT 3*R<1 ) ,V(I J < I > 

0.T600 


F0RMATC1P3E12* 4) 

m 6:5 5 


IFCRTFL. EO»**N0”) G) TO 10 

5 5 60 6 

. 

WR I T F. < "F I L El "> A > N T 

0 3 60 7 

4 

FORMAT Cl 55 

•30 6 1 0 


IvRITFC’FILEl *V5> ( R < 1 0 * VC I) » 1= 1 , NT ) 

03 6 1 1 

5 

FORMAT (FI 6- 8* El 1 • 3 > 

0:3 630 


GO TO, 10 

30 630 


.END 




Reference 


1. J. V. Uspensky, Theory of Equations , p. 175, McGraw-Hill, 
New York (1948). 



APPENDIX D 


MISCELLANEOUS SYSTEM PHOTOGRAPHS 




Figure D-l Receiver and transmitter mounted on azimuth and elevation positioner 
•with dust covers removed. T)nnnAnTw<ni 

REPR0DUC1BIL1 i i w . 
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Figure D-2 Receiver and transmitter mounted on azimuth and elevation 

positioner with dust covers secured. REPRODUCIBILITY OF THE 

ORIGINAL PAGE IS POOR 
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Figure D-3 Laser system power supply and remote 

stations for control of laser and Pockels cell. 









Figure D-4 Electronics console containing Photomultiplier 
tube power supplies and laser power meter 
(on top of console * 
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Figure D-5 Tektronix model 556 Dual Beam Oscilloscope and 
Tektronix C12 Camera. 
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